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lated in chiral perturbation theory at next-to-leading order. We show that the baryon electric dipole 
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1 Introduction 

The neutron electric dipole moment (EDM) is a sensitive probe of CP violation in the Standard 
Model (SM) and beyond. The current experimental limit \d n \ < 2.9 • 10~ 26 ecm [1] is still orders 
of magnitude larger than the SM prediction due to weak interactions. Furthermore, in quantum 
chromodynamics (QCD) the breaking of the U(1)a anomaly allows for strong CP violation, which 
is parameterized through the vacuum angle 9q. Therefore, an upper bound on d n allows to constrain 
the magnitude of 9q. Furthermore, such electric dipole moments are very sensitive to physics 
beyond the SM, see e.g. [2]. Many extensions of the SM in fact lead to larger EDMs than the tiny 
SM predictions, so that any limit of d n leads to bounds on the scale of the new physics. In this 
paper, we concentrate on the CP violation generated by the #-term of QCD. New and on-going 
experiments with ultracold neutrons strive to improve the aforementioned bounds even further, see 
e.g. [3] for a recent review. Furthermore, there are new experimental proposals to measure the 
EDM of the proton and the deuteron in storage rings, which in principle allow for an even higher 
sensitivity than obtained with the instable neutron, see e.g. refs. [4-7]. 
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Besides these challenging experimental activities, first full lattice QCD calculations of the 
neutron and the proton electric dipole moment are becoming available. There exist three different 
methods of calculating the nucleon EDM on the lattice. The EDM can be related to the energy 
difference of the nucleon with different spin alignments in the presence of an external electric 
field [8-1 1]. It can also be obtained by calculating the electric dipole form factor (EDFF) at finite 
momentum transfer q 2 , see the definition in eq. (3.1), and extrapolating to the point with q 2 = 
0[12, 13], 

(i.D 

In addition, the nucleon EDM can also be calculated by analytically continuing 9q to a purely 
imaginary quantity [14, 15], as the QCD action in the presence of the #-term becomes real in 
Euclidean space. For a brief review of these methods, see ref. [16]. We do not discuss here another 
method [17] that relates certain Fock state components of the EDM to the Fock state expansion of 
the magnetic moment (derived in light- front QCD) [18] since it is not clear what these relations 
imply for the observable quantities [19, 20]. 

These lattice studies require a careful study of the quark mass dependence of the nucleon 
EDM to connect to the physical light quark masses. In addition, CP-violating atomic effects can be 
sensitive to the nuclear Schiff moment, which receives a contribution from the radius of the nucleon 
electric dipole form factor (EDFF), see e.g. [21]. It is thus of paramount interest to improve the 
existing calculations of these fundamental quantities in the framework of chiral perturbation theory 
(CHPT). In [22], the electric dipole moments of the neutron and the A were calculated within 
the framework of U(3)lxU(3)j? heavy-baryon chiral perturbation theory and an estimate for 9q 
was given (for earlier works utilizing chiral Lagrangians, see [23-25]). In [26], the electric dipole 
form factor of the nucleon was analyzed to leading one-loop accuracy in chiral SU(2). In that 
calculation, the form factor originates entirely from the pion cloud. The strength of the form factor 
was shown to be proportional to a non-derivative and CP- violating pion-nucleon coupling g-xNN 
that was estimated from dimensional analysis in [26]. In ref. [27], the results of refs. [22, 26] 
were extended to higher order based on a covariant version of U(3)i,xU(3)r baryon CHPT, and 
an expression for g^NN was given in terms of the measurable quantities. For other recent work on 
these issues, see [20, 28]. Here, we extend the studies of ref. [27] to the baryon octet and make 
contact to recent lattice QCD studies. 

Furthermore, the leading contributions to the neutron EDM at finite volume and in partially- 
quenched calculations were considered in [29], and in [30] the leading order extrapolation formula 
using a mixed action chiral Lagrangian is given. Here, we work out the finite volume expression 
for the whole baryon octet at next-to-leading order (NLO). In particular, we show that the for- 
mally suppressed NLO corrections are sizeable, and even dominate for some baryons due to large 
cancellations at leading order. 

The manuscript is organized as follows. The underlying effective chiral U(3)lxU(3)# La- 
grangian is given in section 2. In section 3 we work out the baryon EDFFs and EDMs in the 
infinite volume for varying quark masses. We derive new relations for all baryon EDMs and show 
that these only depend on two combinations of unknown low-energy constants (LECs). Using re- 
cent lattice data, we can give predictions for all baryon EDMs. Then, in section 4, we calculate the 
finite volume corrections for the baryons at NLO and show that these NLO corrections are substan- 
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tial and must be included in any extraction from lattice data. We end with a summary and outlook 
in section 5. Various technical aspects of our calculations are displayed in the appendices. 

2 Effective Lagrangian for strong CP violation 

The most general gauge-invariant and renormalizable Lagrangian for QCD is 

Cqcd = -\G%G a ^ + q(iip-M)q + O^G^G^ (a = 1, . . . , 8) , (2.1) 

with G" the gluon field strength tensor and G® = e flu \ a G a ' Xcr /2 its dual, q collects the quark 
fields of various flavors, is the gauge-covariant derivative, and M. is the quark mass matrix. 
The last term is the so-called #-term, which breaks the P and CP symmetries. It is a consequence 
of the 11(1)^4 anomaly. Because the theta-term is related to chiral U(l) transformations of the quark 
fields (see, for instance, ref. [31]), only the combination 

6 = 8 + &rgdetM (2.2) 

is a measurable quantity. In CHPT, one may treat the #-term using an external field 9{x), and the 
QCD Green functions can be obtained by expanding the generating functional around 9(x) = 9q 
with real quark masses. Under an axial U(l) transformation, one has 

9(x) -> 9(x) - 2N f a, (2.3) 

where Nf is the number of flavors, and a = (an — aL)/2. In the limit of infinitely large number 
of colors N c , the U(1)a anomaly is absent. In this case, the spontaneous chiral symmetry breaking 
of U(3)lxU(3)^, which is a symmetry of the QCD Lagrangian, into U(3)y gives nine Goldstone 
bosons (the SU(3) flavor octet {tt^, it , K^, K°, K°, 7]$} and the flavor singlet r/o). Collecting 
these fields in U(x), this transforms under the axial U(l) transformation as 

U(x) -> e iaR U(x)e- iaL . (2.4) 

Thus, the combination 

9{x) = 9(x) -ilndetU(x), (2.5) 

is invariant under chiral transformations. 

With 9(x), the most general chiral effective Lagrangian which is invariant under U(3) l x U(3)# 
can be constructed. The original construction at order O (5 2 ) can be found in [32], where 0(5) = 
0(M,p, k) with and k denoting the Goldstone boson masses and a small momentum, respec- 
tively. Further, we count 1/N C as O (<5 2 )- Here, we adopt the notation used in [22], whose formu- 
lation is partially based on refs. [33, 34]. The most general chiral effective Lagrangian for mesons 
to second chiral order, complying with the U(3)l xU(3)r symmetry, reads 

C = -V + Vx TrfV^V^] + V 2 Tr$U + + iV s Tr[x^U - xll ] ] 

+V 4 Tr [UVpifi] Tr [uW^U] + V 5 Tr [V^9V^9] , (2.6) 

where x = 2B§M q with M q = diag(m u , m^, m s ) the real quark mass matrix, V M ?7 = d^U — 
iffjU + iUlfj,, and the V^'s are functions of 9(x). 
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In order to use the above Lagrangian, we need to determine the vacuum. Denoting the vacuum 
expectation value of U by Uq, we can decompose U as 



U 



f UoUWU , 



(2.7) 



where 



with 



7T 



V 



7T 

K~ 



(2.8) 



(2.9) 



When #o = 0, the vacuum is trivial, Uq = 1. When the #o angle is switched on, the vacuum is 
shifted, and Uq has to be determined by minimizing the zero modes of the Lagrangian. One may 
parametrize the vacuum as 

U Q = diag (e-^«, e~ iVd , e~^ s ) . (2.10) 
After the vacuum alignment, one obtains the effective Lagrangian [22] 1 



-Vq + Vi Tr 



+ (V 2 + BV 3 ) Tr 



17 + 17+ 



iAV 2 Tr 



+.4V3 Tr 



U + tf 



+ V^Tr 



f/v„,rjt 



Tr 



(2.11) 



where x = 2£?odiag(m u cos ip u ,nid cos ipd, m s cos </? s ), and .A, B are complicated functions of the 
Vi's, see e.g. ref. [22]. Introducing the notation 9q = 9q — J2 q tPq> one ma y expand the V^'s in terms 
of 9q. All of them except for V3, which is an odd function, are even functions of 9q. In the leading 
approximation, A and B are given by 



A 



(2) 



V, 



(0) 



8 + O (<5 4 ) , B 



V* 



(i) 



V, 



(0) 



#0 + O (£ 6 ) • 



(2.12) 



Since 



\/6 

In det U = Tr In f7 = i — 770 , 



the Vi are now functions of 9q + y/Grjo/Fo. The correct normalization of the kinetic terms of the 
Goldstone boson fields can be obtained by requiring 



V 1 (0) = V 2 (0) 



F 2 

_ IT 

T' 



Vk(0) 



12 



W " Ft) . 



(2.13) 



Finally, the quantity 9 can be expressed in terms of the measurable quantity #0 [27], 



70 



(2) 



AM 2 K - Ml 



F* M|(2Mf-M2) 



-1 



70- 



(2.14) 



'in ref. [22], there is one more term i(Vs — BV2)Tr[x (U — U^)]. However, this term vanishes exactly because B is 
defined as V3/V2. 
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One sees that 6$ = O (5 2 ), because Vq 2 ^ is of the zeroth chiral order. 

Similarly, one can also construct the most general effective Lagrangian in U(3)lxU(3)/j 
CHPT for the baryon octet (for a general discussion of effective Lagrangians for the 6>-term, see 
[20]) 

£+ p \ 



B 



V 



1 yo , La 

CO 



(2.15) 



The Lagrangian up to the second chiral order is (only the terms relevant to our calculation are 
displayed; for details, see ref. [22]) 

C+b = iTv[B^[D^B]} -mTr[BB] - jTr[B^ l5 {u^, B}] - ^Tr[B^ Fj [u^ B}} 

+ ^Tr[B^ 5 B}Tr[ Ufl } + b D Tr[B{x + ,B}}+b F Tr[B [x+, B]] + b Tr[BB] Tr [x+] 



/fi 

+44 w' 10 ^- Vo Tr[BB} + i ( w' 13 8 
V6 



W13 



V6 
F 



Vo )Tr[Ba^ 5 {F+,B}] 



(2-16) 



+i(w' u 9 + ^w^r^ojTr [Ba^-ft [F+, B]] , 

where \+ = X+ ~ iA{U — W), and we use the same notation w' l0 = w\q + 2>w\2/2 as in 
ref. [27] with the low-energy constants (LECs) w\q and w\2 defined in ref. [22]. We remark here 
that although there seems to be quite a number of unknown LECs, i.e. wo,w' 10 , wis, w' 13 , wu and 
w[ 4 , only two combinations of these will finally appear in the expressions of the baryon EDFFs. 

3 Baryon electric dipole form factors in the infinite volume 



The electromagnetic form factors of a baryon are defined by 



(B{p')\J^\B{p)) = u(p>) 



1 V F X (q 2 ) 



2 v iF 2 (q 2 ^ 



2ms 



+i (7V75 - 2m B q v ~ls) F A (q 2 ) 



Fs(q 2 
2niB 



u(p), (3.1) 



where F\(q 2 ) and F2(q 2 ) are the P- and CP-conserving Dirac and Pauli form factors, respectively, 
while FA{q 2 ) is the P- violating anapole form factor, and F 3 (q 2 ) if the electric dipole form factor 
(EDFF) which breaks P and CP symmetries. Here, J^ m is the electromagnetic current, the baryon 
mass is denoted by rag, and = p'^ — p^ is the four-momentum transfer. 

In what follows, we will consider the dipole form factor F 3 {q 2 ). The electric dipole moment 
of a baryon is defined as the electric dipole form factor at q 2 = 



d B 



Fa,B(Q) 
2m B 



(3.2) 



Note that similar to the case of the neutron electric form factor, we do not include the normalization 
of the form factor at q 2 = in this definition. 
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(a) 



(b) 



® * 

(c) 




Figure 1. Feynman diagrams contributing to the baryon EDFFs up to next-to-leading order, where (g) denotes 
a CP violating vertex, black dots represent the second order mesonic vertices, filled squares and diamonds 
are the first and second order baryonic vertices, respectively. 
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Table 1. Tree-level contribution to the EDFFs of the octet baryons. 



3.1 Baryon electric dipole form factors up to NLO 

Consider first the tree level contributions to the baryon EDFFs, figure l(a,b). In particular, these 
feature the counterterms w' 13 14 and ^13^4 at O (5 2 ) [22, 27]. The expressions of the form fac- 
tor F^{q 2 )/{2m) of the ground state octet baryons from the tree-level diagrams are collected in 
table 1, 2 where a = lAAV^V^ / (FqF^M^ 2 . One notices that the tree-level contributions to 
various baryons depend only on two combinations of the LECs: awyz + w' 13 and aw^ + w' lA . The 
dependence on the latter may be written in the general form —AQseOo^awu + w' l4 ), with Qb the 
baryon electric charge. 

There are four loop graphs contributing to the baryon EDFFs up to NLO, see figure l(c-f). 
Formally, there are more diagrams. However, due to cancellations they only contribute starting at 

2 The overall sign of the tree-level expression for the neutron EDM given in ref. [27] should be positive. 
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Table 2. Possible loops contributing to the EDFFs of the nucleons and the corresponding coefficients C cc j 
and C c f . The intermediate states of the loops are listed in the second column for each nucleon. 



next-to-next-to-leading order. For details, see refs. [27, 35]. The expression for the sum of loops of 
a given pair of any meson-baryon intermediate state is 

F^tf) _ 8e9 V^ f„ _ [ r ____,^ , (^,_ ^ , „ /f2 f\ T fJ2 _ 2 



a 



cd 



Jmm(q ) + ( 2m(m - rh) + M - — ) JuMmiq , m ) 



2m F% 

+ (C cd + C c{ ) J M m(m 2 )| +0(5 A ), (3.3) 

where m is the mass of the external baryon, M and rh are the masses of the meson and baryon in 
the loops, respectively. We have made use of the fact that the baryon mass difference is O (<5 2 ) , see 
appendix A. Keeping only the leading order, the loop expression does not depend on the baryon 
mass, and reads 

^--•^^W). (3-4) 

The expressions for the loop functions involved in eqs. (3.3) and (3.4), making use of infrared 
regularization [36], are collected in appendix B. 

Different baryons receive contributions from loops with different intermediate states. For in- 
stance, the loops for the neutron can be {7r _ ,p} and {K + , while they are {7r + , n}, {7r°(r/8, i]o),p}, 
{K + , S°(A)} and {K , S + } for the proton. A list of the possible intermediate states and the cor- 
responding coefficients C c d and C e f for the nucleons are given in table 2, where we have defined 

P = {2D- 3w a ) (2b D + 36 + 6u/ 10 ) , (3.5) 

for brevity. The loops and coefficients for the S, A and S hyperons are collected in table 3. The 
loops are divergent. Up to O (<5 3 ) , the divergences can be absorbed into the renormalization of w' 13 



-7- 



and w[ 4 , 



24K (2) 



oy(2) 

w'u = v&(p) + ( 5Db D + $Fb F ) L, (3.6) 

where the divergence is contained in 

L =Cp{^^ [Kto)+r(1,+1] }' (3 - 7) 

with d the number of space-time dimension, and w[ r 3 14 (/x) are the finite parts of w' 13 14 . Notice 
that w' 13 14 (^) depend on the renormalization scale fi through L. The scale dependence cancels 
with that of the loops, and as a result, the final expressions of the EDFFs are scale-independent. 
With the coefficients collected in tables 2 and 3, one can easily obtain the explicit expressions of 
the baryon EDFFs up to NLO. They are listed in appendix C. 

So far, there is no constraint on the LECs w' 13 14 (/x) and 1013,14, which appear at the LO of 
the baryon EDFFs, except for the N c scaling and the naturalness requirement of the effective field 
theory. As mentioned in section 3.1, the tree-level expression of the baryon EDFFs depend only on 
two combinations 

w a {n) = awiz + WiKfj,) (3.8) 

and awu+WiKfj). In principle, these two combinations may be extracted from lattice calculations 
of the baryon EDMs. On the lattice, the calculations can be performed at different quark masses, 
or equivalently pion masses. From a fitting to the pion mass dependence of some of the baryon 
EDMs, one may extract the unknown LECs, and then make a prediction of the other baryon EDMs. 
Furthermore, up to the order O (<5 3 ), 

V (2) (3 

w b (fi) = 3[aw 14 + w[ r 4 (n)} + 2 ° 2 M m (3.9) 

n^STve 

always appear together in the expressions of the baryon EDFFs because of SU(3) flavor symmetry. 
The second term is due to the loops involving the r/o> and its pion mass dependence starts from 
O (<5 4 ) . Here, all the baryon masses in the second term have been replaced by the average mass 
of the baryons m avc since the difference is a higher order effect. Consequently, NLO one-loop 
expressions for all baryon EDFFs are given in terms of just two unknown LECs. 

3.2 Numerical results for the loop contributions of the baryon EDMs 

Since the values of the LECs w' 13 14 (//) and 1013,14 are not known, we will only focus on the 
contributions from the loops in this section, and discuss several relations which are free of these 
parameters in section 3.4. In order to get the numerical results, we use D = 0.804 and F = 0.463. 
From fitting to the baryon mass differences at O (S 2 ), see appendix A, we get bp = 0.068 GeV -1 
and bp = —0.209 GeV -1 . There are other determinations of bo and bp from higher order analysis 
of various baryonic properties, see, for instance, refs. [37, 38]. The difference reflects higher 
order effects, so that they will not be used here. In the large iV c -hmit, one has Fq = F w , and 
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Table 3. Possible loops contributing to the EDFFs of the hyperons and the corresponding coefficients C c d 
and C c f ■ The intermediate states of the loops are listed in the second column for each hyperon. 
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Figure 2. Loop contribution to the neutron EDFF up to LO (left) and NLO (right). The solid, dotted, dashed 
and dot-dashed lines are for the pion mass 138 MeV (physical value), 200 MeV, 300 MeV and 400 MeV, 
respectively. 
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Figure 3. Loop contributions to the EDMs of the neutral baryons as a function of the pion mass. The bands, 
reflecting uncertainties by varying the renormalization scale between /i = M„ and m=;, between solid and 
dashed boundaries are the NLO and the LO results, respectively. The filled circle and square with error bars 
are the lattice data from refs. [11] and [16], respectively. 



we take 92.2 MeV [39] for its value. From an analysis of the 77 — rf mixing in the framework of 
U(3) L xU(3) fi CHPT, it was found that V {2) = -5 x 10" 4 GeV 4 and V 3 (1) = 3.5 x 10" 4 GeV 2 [40]. 
With these values, the neutron EDFF calculated at fj, = 1 GeV is plotted in figure 2, where only 
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the loop contribution is taken into account. One sees sizable NLO effects, especially for higher 
pion masses. One notices that the dependence on q 2 can be well approximated by a linear function, 
which means that a linear extrapolation from finite to vanishing q 2 can be used on the lattice. The 
same is true for the other baryons, as none of them shows a strong g 2 -dependence. 

At the physical pion mass, we get the loop contributions to the baryon EDMs in units of 
KT 16 e6> cm, 

d l °°P = -3.1 ±0.8, 4° OP= 5.6 ±1.0 -6.1 (0 -Gey), 
d l ° op = -2.6 ± 0.4, 4?° p = 3.8 ±1.0 -4.8 (J3 -GeV), 
d l °° p = 0.8 ±0.4, d l °° p = -2.1 ± 0.2 + 4.8 GeV), 
d^ op = -3.6 ± 0.8, d l °° p = -3.7 ± 0.2 + 4.3 (/3 • GeV), (3.10) 

where the uncertainties are estimated by varying the scale \i between the masses of the p and 5. If 
we replace the baryon masses in the /3-term in the baryon EDFF expressions, eqs. (C.2), (C.4), (C.5) 
and (C.8), by the averaged baryon mass m ave = 1151 MeV as that in eq. (3.9), this term contributes 
-5.0 (/3 • GeV) to the proton and £+ and 5.0 (/3 • GeV) to the XT and S~. The difference reflects 
part of the higher order uncertainties. In the following, we will use the averaged baryon mass for 
the /3-term and keep different masses for the other terms. 

In order to compare with the results from lattice simulations, we should study the pion mass 
dependence of the pertinent quantities. We take the physical values for F^, M m . no and the baryon 
masses, since their M^-dependent effects contribute at higher orders. For the kaon mass, we use 

M 2 K = M 2 K + A ^, (3.11) 

where Mk = 484 MeV is the kaon mass in the SU(2) chiral limit with vanishing up and down 
quark masses. The pion mass dependence of the EDMs of the neutral baryons is shown in figure 3, 
where the shaded bands between the solid and dashed boundaries are the NLO and LO results, 
respectively. It is clear that all the EDMs vanishes in the chiral limit. The physical reason is that 
the #-term can be rotated away if any of the quarks are massless, see, for instance, ref. [31], which 
ensures a vanishing value for the #-term induced EDM. 

For comparison, the lattice data for the neutron at M n = 530 MeV calculated in refs. [11] are 
shown, which are calculated on 24 3 x 48 lattice with lattice spacing a « 0.11 fm. We also show 
the lattice data reported very recently in ref. [16] whose uncertainty is smaller. Both the LO and 
NLO loop results agree with the lattice data. 

One notices that the NLO corrections for all the neutral hyperons are dramatic. For the A, 
from table 2 and eq. (3.3), one finds that the loops involving ir + and E~ cancel with those of 7r~ 
and S + exactly. Thus, the remaining LO contributions come from loops involving a kaon. They 
are small as can be seen from 

M 2 Ml 
1 + ln— f = -2.96, 1 + In — f = -0.40 , 

with [i = 1 GeV. Therefore, the LO result for the A EDM is close to zero (as already pointed out 
in ref. [22]). The same happens for the S°. 
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M n [GeV] M n [GeV] 

Figure 4. Loop contributions to the EDMs of the charged baryons as a function of the pion mass. The bands, 
reflecting uncertainties of the LO contributions by varying the renormalization scale between /i = m p and 
ms. The filled circle and square with error bars are the lattice data from refs. [11] and [16], respectively. 
The solid, dotted and dot-dashed lines represent the NLO results evaluated at fi = 1 GeV with /3 = 0,1 and 
— 1 GeV" 1 , respectively. 



The results for the charged baryons are shown in figure 4, where the bands are the LO loop 
results. We choose three different values for the unknown combination j3 to illustrate the NLO 
effects. The results are quite sensitive to the numerical value of /3, especially for the H~ and £ _ . 

3.3 A first determination of the LECs from lattice results 

We may use the neutron and the proton EDMs at unphysical quark masses to determine the two 
combination of parameters w a (fi) and w^^fi). Using the lattice data at M w = 530 MeV [16] 3 for 
the neutron EDM, one may determine the counterterm combination w a {n) 

w a (l GeV) = (-0.01 ± 0.02) GeV" 1 , (3.12) 

where the uncertainty merely reflects the uncertainty in the lattice calculations. Using the data for 
the proton, we get 

w b (l GeV) = (-0.40 ± 0.05) GeV" 1 . (3.13) 

3 As this pion mass value is at the edge of the range of applicability of our approach, lattice data at lower pion mass 
are urgently called for for a more reliable determination of the LECs. 
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Baryon 


p n £+ S° S- A ~° S" 


Combination 


-{w a + w b ) 2w a -(w a + w b ) -w a w b -w a w a 2w a w b -w a 



Table 4. Combinations of the unknown parameters appearing in the EDFFs of the baryons up to NLO. 



With these determinations, the EDMs at the physical pion mass can be predicted including both 
the tree and the loop contributions. The ones for the neutron and proton in units of 10~ 16 e 6$ cm 
are 

d n = -2.9 ±0.4 ±0.8, d p = 1.1 ± 0.5 ± 1.0. (3.14) 

Here, the first uncertainty reflects the uncertainty in the determination of w a (l GeV) and w b (l GeV), 
and the second one corresponds to varying the scale fi between the p-meson mass and m=. Com- 
paring with the loop results given in eq. (3.10), one sees that the neutron EDM at the physical pion 
mass is dominated by the loops. For the proton, the loop contributions have similar size as, albeit 
slightly larger than, the tree level terms. Combining these two errors, we have d n = —2.9 ± 0.9 
and dp = d- p = 1.1 ± 1.1 (in canonical units). Using the experimental upper limit of the neutron 
EDM, \d n \ < 2.9 x 10~ 26 e cm [1], the 6q angle is constrained to be 

|0 O | < 1-5 x 10- 10 . (3.15) 

Similarly, we can predict the hyperon EDMs (again in units of 10~ 16 e 6>o cm) 

d A = -2.5 ± 0.2 ± 0.4, d 2+ = -0.7 ± 0.5 ± 1.0, 
d s o = 0.7 ±0.2 ±0.4, d E _ = 2.2 ±0.5 ±0.2 
d s o = -3.4±0.4±0.8, d B - = 0.6 ± 0.5 ± 0.2. (3.16) 

3.4 Counterterm-free relations 

In this section, we will derive relations that are free of the unknown LECs up to NLO. These 
relations can e.g. serve as checks of lattice simulations for varying quark masses. In fact, due the 
SU(3) flavor symmetry, there exists an exact relation among the EDFFs of the £ hyperons. From 
table 2, one deduces that the loops involving a pion or an ry 8 ( ) for the S + cancel with those of 
the As mentioned before, the pionic loops of the S° cancel with each other. Furthermore, 
there are also cancellations for the kaonic loops among the £ hyperons. These cancellations occur 
when isospin is a good symmetry. Therefore, in the isospin symmetric case, which is considered 
throughout, one has 

^3,e+ + ^3,e- - 2F 3iE o = O (<5 4 ) . (3.17) 

This relation is exact up to NLO including the tree-level contributions. Similarly, there is another 
relation which is exact only at LO 

F 3)E o+F 3) a = 0(<5 3 ), (3.18) 

since the LO loop contributions are independent of the baryon mass, and they cancel each other in 
the sum. Moreover, the above combination is also independent of /3 which combines the unknown 
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Figure 5. Counterterm-free combinations of the baryon EDMs as a function of the pion mass. 



parameters wo and w' 10 , so that it can be calculated parameter-free. Taking q 2 = 0, the O (<5 3 ) 
expression is rather simple 

(2) □ 

dso +d A = - f*/ 4 J° {M 2 K - M 2 ) (Fb 2 D + 2Db D b F + 3Fb 2 F ) + O (5 4 ) . (3.19) 

In fact, there are more combinations free of the counterterms. As mentioned before, up to 
NLO, there are only two combinations of unknown constants, w a {(j) and Wb{p), in the expressions 
of the EDFFs of the octet baryons defined in eqs. (3.8) and (3.9). It is easy to find combinations 
free of any unknown parameters up to O (<5 3 ) utilizing the results collected in table 4. In addition 
to the two in eqs. (3.17) and (3.19), they include 

F3,p-F 3 ^+, F 3:n -F 3jS o, F 3>E --F 3 p-. (3.20) 

Because the counterterms have been cancelled out, these combinations are also finite and hence 
independent of the scale //. However, different from eq. (3.19), their LO loop contributions do not 
vanish. The expression for the F 3;H — F 3 =o at q 2 = reads 



d n — dwo 



eV o 2) ^\ mh ,( , M 2 K M^-M K \ 

(Db D + Fb F ) 2 In — + vr 1 



+^{M 2 K -M 2 ) (Db 2 D + 2Fb D b F + Db 2 F ) +0(5 4 ). (3.21) 
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The relations for the charged baryons are more complicated, and are not shown here (they can 
be obtained from the expressions given in appendix C). The pion mass dependence of the four 
combinations in eqs. (3.18) and (3.20) are shown in figure 5. Their values at the physical pion mass 
in units of 10~ 16 e 9q cm are 

d s o + d\ = — 1.8, d- p — <i s + =1.8, 

d n - dwo = 0.5, d E _ - d s - = 1.6 . (3.22) 

Note that although the first combination starts from one order higher than the others, their numerical 
values are of similar order of magnitude. 

4 Finite- volume corrections in the p-regime 

On the lattice, calculations are performed in a finite volume. As a result, the continuum mo- 
mentum spectrum becomes quantized. Taking periodic boundary condition for all three spatial 
dimensions forms a torus. If the volume is L 3 , the momentum takes values of 2nn/L, with n a 
three-dimensional vector of integers. One easily sees that any integral over the spatial components 
of momentum in the infinite volume should be replaced by a summation for the momentum modes. 
For instance, the two-point scalar loop integral needs to be changed as follows, 



d 4 k 1 i ^ f dk° 



(2vr) 4 (k 2 - mf)[(k + q) 2 - m 2 ] L 3 ^ J 2vr (k 2 - m 2 )[(k + q) 2 - m 2 } 



(4.1) 



The finite volume corrections to a quantity Q is defined as the difference of Q evaluated in a 
finite and infinite volume, 

S L [Q] = Q{L) - Q(oc). (4.2) 

As we already mentioned, at chiral limit the baryon EDMs should vanish. However, the current 
lattice data at rather large pion masses do not show a decreasing behavior yet [10, 1 1]. Thus, results 
at smaller pion masses are necessary. However, in the section, we will show that the finite volume 
corrections at small pion masses to most of the baryon EDMs are quite large, and it is necessary to 
include the NLO contributions. The decomposition of the matrix element of the electromagnetic 
current in the form of eq. (3.1) assumes Lorentz invariance. However, on discretized lattice, Lorentz 
invariance in infinite volume is lost. Furthermore, the external momentum is also quantized in 
a finite volume. These effects should be taken into account when calculating the finite volume 
corrections to the form factors [41, 42]. Since we are not aiming at a very precise calculation of 
the finite volume effects, we will simply assume Lorentz invariance in the following. 

As shown in ref. [43], one can use the same LECs as in the infinite volume, and the finite 
volume corrections come only from loops. This may be understood easily, since the finite volume 
effects are long-distance physics, while the LECs reflect the short-distance physics. For the same 
reason, the finite volume corrections should not depend on the choice of method to regularize the 
ultraviolet divergence. This means that one may change the upper bound of the Feynman parameter 
integration in eqs. (D.5) and (D.7) in the infrared regularization from oo to 1 as already noticed in, 
for instance, ref. [44]. In the p-regime, where M W L S> 1 [43, 45], the finite volume corrections are 
exponentially suppressed. 
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Figure 6. Left: The ratios of the finite volume corrections to the loop contributions in the infinite volume 
for the neutron EDM as a function of L. Right: The ratios of the finite volume corrections at NLO (LO not 
included) to those at LO. The solid, dotted, dashed and dot-dashed lines are for the pion mass being physical, 
200 MeV, 300 MeV and 400 MeV, respectively. In the left panel, for each value of the pion mass, two lines 
are plotted with the upper and lower ones representing the NLO and LO results, respectively. 



In ref. [29], finite volume correction to the neutron EDM was calculated at LO in SU(2) heavy 
baryon CHPT. The resulting expression is given by [29] 

x r^LOi 9Aae9 m u m d ^ 

6 L [d n ]= — — ^ ^jr (L^|n|), (4.3) 

71^0 

where a is a coefficient of the CP-violating NNtt vertex. 4 From eqs. (C.l) and (D.4), our result 
for the finite volume correction to the neutron EDM at LO is 

^K°] = Y,[(D + F) (b D + b F )K (LM*\ft\)-(D-F) (b D - b F ) K {LM K \n\)\ , 



where we have used [22] 



(4.4) 

- F 2 M 2 
8V {2) 

which may be obtained from eq. (2.14) considering M w <C Mk- Comparing the SU(2) part with 
eq. (4.3), we can identify 

a = -^{bo + b F ). 

The left panel of figure 6 shows the ratios of the finite volume corrections to the loop contri- 
butions in the infinite volume for the neutron EDM as a function of L. The results are obtained 
with fj, = 1 GeV. In order to show the impact of the NLO corrections more clearly, we plot the 
ratios of the finite volume corrections at NLO (LO not included) to the LO in the right panel of 
figure 6. One sees that at the physical pion mass the NLO contributions increase the LO finite 
volume corrections by about 20%. This correction is much larger for a pion mass of 400 MeV. At 
this point, one might worry about the convergence. In fact, the neutron EDM in the infinite volume 



4 a is the a used in ref. [29]. 
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L [fm] L [fm] 

Figure 7. Left: The ratios of the finite volume corrections to the loop contributions in the infinite volume 
for the EDMs of neutral hyperons as a function of L. Right: The ratios of the finite volume corrections of 
NLO to those of LO. The solid, dotted, dashed and dot-dashed lines are for the pion mass being 138 Nev 
(physical value), 200 MeV, 300 MeV and 400 MeV, respectively. 

has a nice convergence property, as can be seen from the left panel of figure 3. Moreover, one 
can show that the ratio tends to increase with increasing L. Although both the LO and NLO finite 
volume corrections decrease exponentially, the LO one decreases faster. Using eqs. (E.2) and (E.4) 
from appendix E, we get in the limit L — > oo, 




(4.6) 



The finite volume corrections to the neutral hyperons are shown in figure 7. It is obvious that the 




L [fm] L [fm] 



Figure 8. The ratio of the finite volume corrections to the loop contributions in the infinite volume for the 
proton EDM as a function of L. The first graph presents the LO results, and the others are the NLO results 
for different choices of (3. The solid, dotted, dashed and dot-dashed lines are for the pion mass of 138 MeV, 
200 MeV, 300 MeV and 400 MeV, respectively. 



effects are very small for the A and TP. The reason is that the pionic loops with positively and 
negatively charged pions cancel each other exactly for these two baryons. This also makes the 
NLO corrections much larger than the LO ones at fj, = 1 GeV, consistent with the infinite volume 
results shown in figure 3. The situation for the H° is similar to case of the neutron. 

Although the parameter combination Wb{n) can be determined from the lattice data of the 
proton EDM, (3 is still unknown since it always appear together with the LECs wu and w'/^fi). 
However, the finite volume corrections for the charged baryons depend on (3, but not on W14 and 
w'/^fi). It is therefore possible to extract the value of /3 from the volume effects of the baryon 
EDMs. In figure 8, we show the ratios of the finite volume corrections to the loop contributions 
in infinite volume to the proton EDM for different chosen values of /3. At this point, we want to 
emphasize that the EDMs of all the charged baryons depend on the same combination of LECs, 
i.e. p. Once it is determined from one baryon, it can be used for predicting the finite volume 
corrections to the others. 



5 Summary 

In this paper, we extended the previous calculation [27] of the nucleon EDFFs and EDMs to the 
ground state baryon octet up to NLO in the framework of U(3)l x U(3) r chiral perturbation theory. 
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Our main findings can be summarized as follows: 

1) We have shown that the complete one-loop expressions for the baryon EDFFs and EDMs 
depend on two combinations of unknown LECs only, w a and Wb, cf. eqs.(3.8,3-9). In case of 
the charged baryons, the combination Wb{fi) combines two LECs from the tree graphs and 
one LEC that appears only in loops. 

2) We have shown that the NLO corrections are large for the neutral hyperons A and E°. This is 
due to a suppression of the LO contributions based on exact cancellations between loops of 
positively and negatively charged pions. For the charged baryons, we find a strong sensitivity 
to the LEC combination Wb- 

3) We have derived a set of relations between various EDMs that are free of unknown LECs. 
These can be useful for future lattice simulations of baryon EDMs. 

4) Based on recent lattice results for the neutron and proton EDMs at = 530 MeV, we could 
pin down the two LEC combinations w a and Wb- Based on this, we can predict the baryon 
EDMs at the physical pion mass. In particular, we find d n = — 2.9 ± 0.9 and d p = 1.1 ± 1.1 
in units of 10~ 16 8q e cm. 

5) The finite volume corrections to the baryon EDMs in the p-regime are also studied. Because 
the loops contribute from LO, it is found that the finite volume corrections are huge for all 
the baryon EDMs except for the A and S° for the pion mass close to its physical value. For 
the neutron, the finite volume correction is about 10% at M n L = 4. 

Note that the calculation of the finite volume corrections in our paper assumes Lorentz invari- 
ance, and additional subtleties in the presence of external electromagnetic fields, see [46], are not 
taken into account. For a precise calculation of the finite volume corrections of the baryon EDMs, 
separate analyses need to be done for the lattice calculations using the form factor method, the 
spectrum method with external field and methods with the twisted boundary conditions. 

More lattice results of the nucleon EDMs are expected to come out soon. Calculations of the 
hyperons are welcome based on our analysis. With these upcoming lattice data, one can determine 
the LECs more reliably, and a next-to-next-to-leading order calculation would be feasible and de- 
sirable in view of the bad convergence for some baryons. Furthermore, it would be interesting to 
see how the LEC-free relations get modified at higher orders. 
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Baryons 


N £ A E 


K 

b K 


b F b D -b D /3 -b F 
b D -b F Ab D /3 b D + b F 



Table 5. Coefficients b^ t K m tne NLO baryon mass formula. 



A Baryon masses up to NLO 

From eq. (2.16), it is easy to obtain the baryon masses up to O (S 2 ), which may be written as 

m B =m- 2b {Ml + 2M 2 K ) - 4 (M<? + b K M 2 K ) + O {5 3 ) , (A. 1) 

where we have expressed the quark masses in terms of the meson masses at the lowest order, and 
b^,K are listed in table 5. 

B Loop integrals in infrared regularization 

Loop integrals in the infinite volume used in this paper are collected in this appendix. Let us denote 
the masses of the light meson and baryon in the loop by M and m, respectively. The momentum 
of the external current is q^, and the momenta of the external baryons are denoted by p^ and 
p'/j- = _|_ qt* t respectively. We further define 

which contains the divergence at the spacetime dimension d = 4. 

The 2-point mesonic-loop function when the two mesons have the same mass reads as 



(27r) d {k 2 -M 2 + ie) [{k + q) 2 - M 2 + ie] 

— ^ I In — s 1 — a In 

16vr 2 \ fi 2 a+1 



2L +T7^o ( ln — -l-o-ln^-r ) , (B.l) 



where a = ^\-^M 2 jq 2 . 

When there is a baryon in the loop, the loop evaluated in the normal dimensional regularization 
spoils power counting [47]. The power counting can be restored in the heavy baryon formalism of 
baryon CHPT [48, 49]. There are also covariant formalisms, see ref. [50] for a recent review. 

Here, the infrared regularization [36, 51] will be used. In the infrared regularization, the inte- 
gral is separated into an infrared regular piece and an infrared singular piece. The infrared regular 
piece can be expanded analytically in the chiral expansion, and absorbed into the counterterms. 
Only the infrared singular part of the loop integrals will be shown in the following. 

The two-point scalar loop integral involving one meson and one baryon is 

d d k 1 



JMm{p 2 ) 



(2Tr) d (k 2 -M 2 + ie) [(p - k) 2 - m 2 + ie] 

M 2 \ r-— n „ — ^ / U 



oj r 1 



p 2 32tt 2 p 2 



ln — g — 1 I + 2\/ X(p 2 ,m 2 ,M 2 ) arccos 
/ 



2pm / 



(B.2) 
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where ui = p 2 — m 2 + M 2 , p = yfp 2 , and X(x, y, z) = x 2 + y 2 + z 2 — 2(xy + yz + xz) is the 
Kallen function. Noting that the external baryons are on shell, p 2 is given by the mass squared of the 
external baryon. Thus, with the baryon mass expressions in appendix A, one has p 2 



m 



and lo = O (<5 2 )- Expanding the arccos function, we get 



JMm(p 2 



UJ 



M 

lQirp p 



L + 



1 



32^ 2 



1 + In 



7? 



+ (5 3 



+ 0(5 2 ) 



(B.3) 



where the first and second term are of 0{5) and O (<5 2 ) , respectively. 

The three-point loop integral with two mesons and one baryon can be written as 

d d k 



(2ir) d (k 2 -M 2 + ie) [(k + q) 2 — M 2 + ie] [(p - k) 2 - m 2 + i e] 



d 



dM 2 



(B.4) 



where M 



f2 



M 2 + x(x - l)q 2 , 



p M + xq^. In the following, p = p will always be 



assumed, then one has p 2 = p 2 + x(x — l)q 2 . The leading chiral order of this integral is 0{5 1 ). 
Keeping terms up to O(5 ), the analytic expression for the infrared singular part is 



JMMm(q 2 ,p 2 



1 



Wrrp^/— q 2 
L I 
p 



arctan 



2M 



+ 



2 + 32vrV 



1 +ln 



M 2 
V 2 



+ 



2oj - q 2 , a-1 



q 2 a 



In 



a + 1 



+ 0(6), (B.5) 



where the first term is of order 0(5^ 1 ), and the terms in the square brackets are of order O(5 ). 

For completeness, we also give the expression for the loop with one meson and two baryons, 
though not used in our calculations. Similar to JMMm(q 2 ,P 2 ), this loop can also be worked out 
from the two-point loop with one meson and one baryon, 

f rl d h 

i\J> ,P 



J Mr 







(2ir) d (k 2 -M 2 + ie) [(p - k) 2 - m 2 + i e] [(p' - k f 



dxJ M m(p 2 



dm 2 jo 

where fh 2 = m 2 + x(x — l)q 2 . The analytic expression up to NLO is 



Jm mm (P ,P' 



1 



2 + 32tt 2 p 2 



1 + ln 



M 5 

7? 



+ 



LO 



6Airp 2 mM 



+ 0(5 2 



m z + i ej 

(B.6) 

) , (B.7) 



where the terms in the square brackets are at LO, and the second term is at NLO. Up to NLO, this 
loop does not depend on the momentum transfer q 2 . 



C Expressions for the baryon EDFFs up to NLO 

Up to NLO, which is O (^ 3 ), the baryon EDFFs contain two parts: the tree-level expressions are 
given in table 1, and the loop contributions are obtained by taking the LO terms of the baryonic 
loops JMMm(q 2 , w 2 ) and JMfh(q 2 )- The explicit expressions of the EDFFs of the baryon octet up 
to NLO are given in the following: 
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• For the neutron, 



2m n 



V£ 2) e9 



* (2Aff - g 2 ) 

H arctan 

2mArV-g 2 

vr (2M 2 K -q 2 



tt 2 F$ 
2M T 



(D - F) (b D - 6 F ) 



1 l M -o. l a *~ X 
1 - In — v + o-jr In — 

1 — In — + ok In 



2m^ 



8 (6 D - hp) (Mi - M 2 ) j arctan 



2M^ 



0-JS- + 1 
(C.l) 



with a = 144y (2) V 3 {1) /(F F 7r M w ) 2 , and a n{K) = yjl - AM 2 (K) /q 2 . Here, we have used the 

SU(3) mass splitting for the baryon masses — mjy = 4 (bo — bp) {M K — M 2 ) + (<5 3 ). 
Similar relations derived from the Lagrangian eq. (2.16) will be used in the following. 



For the proton, 



2totv 



r(2). 



^ [a (1013 + 3^i4) + + 3<(/i)] - \ 6(D + F) (b D + bp) 



M 2 

1 - In — f + In 



<t w - 1 , 3ttM w tt (2M 2 - g 2 ) / 



6tt 2 F4 

2 J2> 



arctan ■ 



+4 (Dbp> + 3Fb F ) ( 1 - In # + *k In + ^ 



47T a/-? 2 
H — ; arctan 



2M X 



+8 {M 2 K -M 2 )[F (b 2 D + 36|) - (6 D - 36 F ) 



+- 



7T 



2F 2 



6(D - F) (b D - bp) M K + (D- 3F) (b D - 3b F ) M m + -^pM; 



(C.2) 



where /3 = (2D - 3w ) (2b D + 3b + 6w' 10 ) 
• For the S°, 

4 



2m s 



T/( 2 ) a ( / 11x2 

-e0 o Hi 3 + < 3 (M)] " ^TpT + F M ( 1 " ln 



, UK — 1 \ TT 

+(Jk ln — r I H ; 7 arctan 



a K + lJ ^f- q 2 



2M K 



Db F + Fb D 
2ms 



;2M 2 - q 2 ) 



+8 (M 2 K - M 2 ) (Fb 2 D + 2Db D b F + Fb 2 F ) 



(C.3) 



-22- 



For the E+ 



3,S+ 



2ms 



--e6» [a («Ji3 + 3u>i 4 ) + w'i^/x) + Su/^/x)] 



^o (2) e^o 
3vr 2 Fi 



2(Db D + 3Fb F ) 



1 - In — + Utj- In + 



ov + 1 m s 



+3(D + F) (6d + 6f) 
vr (2M 2 K -q 2 



l-ln^ + ^ln^+^ 



A* 



cr^r + 1 m s 



+ 



T 3 ? 

2^ 



: arctan 



2ms 



+ 8 (&£> + 6 F ) (Mr- - M^) J arctan 



2-Mtt 



(Db D + 3Fb F ) 
2ms 



2M^ 

32 

(2M| - g 2 ) + y Dft 2 , (Mi - M*) 



+ - 



7T 



m s 



GFbpM-x + 3(D - F) (6 D - 6 F ) M K + 2Db D M m + -^/3M„ 



.. (C.4) 



For the S" 



^3,E-(g 2 ) 

2ms 



-e6» [a (u>i 3 - 3iui 4 ) + io^//) - S^i^A*)] 



3vHi^ ^ y ^ a n + 1 ms / 



2 - arctan-^ 



2M,r 

+3(D - F) (b D - b F ) 



(Dbo + 3Fb F ) 

2ms V * 

1 - In — ^ + ctr- In H 

ok + 1 ms 



, 2 ) + |DJ^ (Ml - M 2 ) 



+ 



7T 



2M| - q 2 
2m s 



+ 8 (6 D - b F ) (Mk - Ml) j arctan 



+- 



7T 



ms 



2M X 
F 2 

1 TT 



6Fb F M w + 3(D + F) (b D + 6 F ) Mr- + 2Db D M m + ^^M w 



(C5) 



• For the A 
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• For the S°, 
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• For the E , 
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The expressions for the EDMs can be easily obtained by noticing 

Bm^.-J, k -4= arctan ^ - 1 



D Finite volume corrections to loops 



2M 2M 



This appendix is dedicated to finite volume corrections to the two- and three-point loop integrals. 
Let us consider the scalar two-point loop 



d A k 
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(D.l) 



where m 2 2 = x(x — l)q 2 + (1 — x)m 2 + x m 2 .. If the loop involves both a meson and a baryon, mi 
should be replaced by the mesonic mass, and the upper bound of the integration over the Feynman 
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parameter x should be replaced by oo in infrared regularization. Performing the contour integral 
over k°, one gets 

1 f 1 , f d 3 k 1 



w« )—iJ^J (S?/ t - s+ ^w {D2} 



Using the formula 

1 r d 3 k \ 1 ^ 2^a 3 - 2 ^ v %- 3 (£a|n|) 

L s2^ J (27r) 3j ^ 2 + fl2 y ^/^(j) |f (La|n|)3/2-i 



n 



derived in ref. [52], we get the finite volume correction to the loop J mi m 2 ((/)> 

1 Z" 1 

^L[Jm im2 (^ 2 )] = -g-2 / dx^2K (Lm 12 \n\), (DA) 

with K v (z) a modified Bessel function of the second kind. 

From eq. (B.4), it is easy to write down the correction to the three -point loop integral, 

,2 „2m _ ° / j_jr_r t_ 



$L [JmMwi (q ,P )] = dx5 L [J Alm (p )] 



y^n 2 [ dx ^ dy^-K^zx), (D.5) 
"777. JO io z l 



ISvr 2 — .0 JO 

where #i = L|n| \y(y — l)p 2 + !/m 2 + (1 — y)M 2 ~\ . We have made use of the integral repre- 
sentation of K n [z) [53] 

^ , , T(n + l/2)(2z) n f°° , cost 

A - w= ^ /"Wr^ (D ' 6) 

which is valid for Ren > — 1/2, | arg z\ < it/ 2. Similarly, one gets 

-2 /•! /-oo 



<5L[^Mmm(p 2 ,p /2 )] = 7^ V]n 2 dx dy—K 1 (z 2 ), (D.7) 
1,177 ^ :a 

where ^2 = L\n\ [y(y — l)p 2 + ym 2 + (1 — y)M 2 



n^O 



E Asymptotic expansion of finite volume corrections 

For simplicity, we will focus on the asymptotic expansion of finite volume corrections with q 2 = 
and p 2 = m 2 . The general case can be treated similarly. For \z\ — > oo, one has 



1 + 



(E.l) 



For L — > oo, we can consider the term with \ ft\ = 1 only in eq. (D.4). There are 6 possibilities, so 
that for L — > oo, 

q f= -LM 
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For the three-point loop integral, 



j2 r^T poo 



,-LC(y) 



[LC{y)\ 



3/2' 



(E.3) 



with C(y) = [y 2 m 2 + (1 — ^M 2 ] 1 / 2 . The leading term of the above integral over the Feynman 
parameter y can be worked out using Laplace's method. For very large L, the integral receives 
contributions mostly from the neighborhood of the minimum of C (y), which is at y c = M 2 / (2m 2 ). 
Thus, 



h[JMMm(0,m 2 
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8tt 2 V 1 
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(E.4) 



The asymptotic expansion for <5l [ Jjumm ("i 2 , m 2 )] is the same. 
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